Particle Production at High Energy and Large Transverse Momentum - "The Hybrid 

Formahsm" Revisited. 



Tolga Altinoluk and Alex Kovner 
Physics Department, University of Connecticut, 2152 Hillside Road, Storrs, CT 06269-3046, USA 

(Dated: January 13, 2013) 

We revisit the "'hybrid formalism"' for particle production used recently to study saturation 
effects in single hadron multiplicities at forward rapidities at RHIC and LHC. We point out that at 
leading twist there is an extra contribution to the formulae used so far, which corresponds to particle 
production via inelastic scattering of the projectile partons on the target fields. This contribution 
' is expected to be small due to kinematics at very forward rapidities/very high transverse momenta, 

^Nj , but should be significant at high momenta and very high energies. This contribution is expected to 

be most affected by saturation effects and is therefore an interesting object of study in the context 
' of possible onset of saturation at high energies. 
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In the last several years a large amount of RHIC data has been interpreted in terms of perturbative saturation 
physics [1], or the so called Color Glass Condensate (CGC)[2]. Although this interpretation is philosophically very 



simple and appealing, alternative interpretations are also available and it is important to be able to differentiate 
f*) |. between them. One would like to understand to what extent the data really unambiguously supports the idea of 
^ ' saturation. One of the problems we are faced with in this regard, is that the saturation based calculations, although 
r^ .. in principle rooted in controlled perturbative approach to fundamental QCD physics, in practical implementations 
rely on phenomenological or semi - phenomenological ansatze and variety of shortcuts. As a result it is sometimes 
difficult to understand what features of the theoretical results are genuine and robust predictions of saturation, 
^ ' and which are model dependent transient features. The saturation based calculational techniques have advanced 
\ considerably during the last two-three years. In particular large part of next to leading corrections [3] is now taken 
i into account in calculating the evolution of gluon density to high energy. This allowed for good fits to the HERA 
DIS data at low x[4]. Still there is a lot of room for improvement, and in our view this improvement is absolutely 
essential in order for comparisons with data to have clear value. It is especially important at this time, since lately 
' many aspects of the RHIC (and LHC) data have been analyzed in the framework of saturation physics[5]. These 
include the single inclusive particle production[6], [7] and the two hadron correlations at forward rapidity [9], [10] in 
dA collisions, and more recently description of single particle spectra[8] and attempts to understand [11] ridge in p-p 
collisions at LHC [12]. Although qualitatively the agreement of saturation based approaches with data is heartening, 
the calculations have still a long way to go before they can achieve the level of reliability comparable to that of, for 
example standard perturbative QCD approach. 
$^ , The present paper is a modest contribution to improvement of one particular calculational approach addressing one 
particular piece of data - particle production in deuteron - gold collisions at forward rapidities. RHIC experiments 
observed strong suppression in the particle production in dA at forward rapidities. The "'state of the art"' saturation 
calculation of this effect appears in [7]. Although the data is described quite well, there are some peculiarities to 
the results of [7]. First, a very small iC- factor is required to fit the overall magnitude of the production of neutral 
pions, while no K - factor is required to fit the charged hadron multiplicity. Secondly and perhaps more worryingly, 
the suppression in the theoretical curves of [7] when extended to LHC energies persists to extremely high transverse 
momenta, where one expects perturbation theory to be long applicable and RdA to be equal to one. 

The calculations of [7] are based on the '"hybrid formalism"' of [13]. In this approach the wave function of the 
projectile at large values of x is calculated perturbatively, without soft approximation, while the scattering of the 
projectile partons on the target fields is treated in the eikonal approximation. The exact treatment of the projectile 
function is of course necessary to describe particle production at forward rapidity, since these partons cannot be in 
any way considered soft. It has another advantage over the soft approximation in the projectile wave function since 
it satisfies energy (longitudinal momentum) conservation for the incoming projectile. The energy conservation in the 
scattering process is still violated of course, since the recoil (and radiation) in the scattering event itself is not taken 
into account in the eikonal approximation. Although energy conservation must be very important at large values of 
xp and its effect has to be understood to make sure the treatment is consistent, we have nothing new to add to this 
point. In this paper we will revisit the derivation of inclusive particle production within the hybrid formalism per se 
relaxing only the collinear approximation made in [13]. Our goal is to identify the terms which where omitted in [13] 
but may nevertheless be important when the transverse momentum of produced particles is significantly higher than 
the saturation momentum of the target. As we will show, such terms, which do not correspond to collinear emission 
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of the incoming projectile partons do indeed exist and contribute at leading twist. These terms have a simple physical 
interpretation and also have a simple form amenable to numerical implementations. It is out hope that including 
these terms in numerical calculations will improve the quality of the results. 

The expression derived in [13] and used in [7] for particle production has a very intuitively appealing and simple 
form 



dN 1 dz 



d^kdr] (27r)2 J^^ 



xifgixi,Q^)NAix2, -,b = 0)Dh/giz, Q) + ^gXifgixi,Q^)Np{x2, -,b = 0)^^/,(z, Q) 



(1) 



where iV^(j')(fc) is the Fourier transform of the forward scattering amplitudes of the adjoint (fundamental) dipole. It 
describes the process whereby incoming low pt partons scatter on the target independently of each other, acquiring 
large momentum kr in the process, and subsequently fragment into observed hadrons. This process is certainly the 
origin of large part of produced particles. 

However there is another physical mechanism which produces large fcr particles in leading twist, whereby high pt 
particles preexisting in the wave function of the incoming projectile scatter with only a small momentum transfer 
from the target. The soft scattering is nevertheless enough to decohere the incoming partons from the rest of the 
wave function so that they materialize as on shell particles in the final state. The high pT partons in the projectile 
wave function arise due to DGLAP splitting of very forward partons. The scattering process is essentially just the 
inelastic scattering of the forward projectile partons with emission of gluons (or quarks/antiquarks)[14]. As we will 
show explicitly within the hybrid formalism, this mechanism of production is equally important as the one taken 
into account in eq.(l) when the saturation momentum of the target is small. When Qs is large, this contribution is 
somewhat suppressed, but may still be quantitatively quite large. Parametrically, while the contribution of eq.(l) 
is roughly proportional to In j^^, the additional inelastic scattering contribution scales like In^. It is thus only 
suppressed for kx ~ Qs when Qg <^ Aqcd, and even then the suppression is merely logarithmic. Given that for 
RHIC data Qs/^qcd ~ 5, it seems prudent to keep this contribution in numerical calculations. 

It is quite clear that taking into account the inelastic mechanism must bring the calculation of particle production 
into agreement with the perturbative result at large pr- Thus we hope that inchiding this contribution will bring RdA 
close to unity at reasonable values oipT- It is also interesting to note that the final states of the inelastic scattering 
are quite different from those of the elastic one. The elastic piece is dominated by quarks in the final state, while the 
inelastic one contains comparable number of quarks and gluons. Since the fragmentation functions of quarks are very 
different, we expect it to affect the relative magnitude of neutral pion and charged hadron production and thus be 
relevant to the problem of a very small K factor for neutral pion production encountered in [7]. Whether it helps or 
makes things worse remains to be seen. Naively one expects gluons to fragment predominantly into neutral mesons, 
and thus the problem of the K factor may become even more acute, since the neutral to charged hadron ratio is likely 
to increase after including the inelastic contribution. 

We also note that it is the inelastic term which is especially sensitive to the saturation effects. The wave function 
of the incoming hadron knows nothing about saturation by itself. The effects of saturation come entirely from the 
distributions of the target. The target fields are directly affected by saturation at momenta k < Qg. The elastic 
scattering probes the large momentum component of target fields, equal to the final momentum of the produced 
parton. Thus as long as pt > Qs, this part of hadron production should be less affected by the saturation effects 
and one could expect that its dependence on energy and atomic number stems from perturbative physics. Any non 
leading twist scaling then presumably comes from effects of a possible " 'nonperturbative" ' initial condition propagated 
to higher momenta via perturbative evolution. The inelastic scattering contribution on the other hand probes the 
target fields at kx ^ pr which includes the region kx < 0{Qs)- It is this region of momenta which is strongly affected 
by target saturation effects. Thus if one neglects the inelastic contribution, one also severely limits ones options of 
studying ciffects of saturation. 

The structure of this paper is the following. In Sec. 2 we revisit the derivation of particle production in the hybrid 
formalism. For simplicity of exposition we discuss in this section the pure glue theory. In Sec. 3 we generalize our 
derivation to full QCD with quarks and antiquarks. We present a simple formiila for the inelastic parton production 
contribution in terms of the gluon and quark distribution functions. We conclude in Sec. 4 with discussion of our 
results. 



II. A PROTOTYPICAL CALCULATION - GLUON PRODUCTION 



In this section we derive the expression for gluon contribution to hadron production in the hybrid formalism. We 
will include the quark and antiquark contributions in the following section. Our approach in the formal sense is similar 
to that of [15], although like in [13] we are not approximating the gluon splitting function by its low x limit. This will 
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give us a possibility to compare our results with the kr factorized formula which arises very simply in the approach 
of [15], to get some intuition from the simple kr factorized expression and also to see the similarities and differences 
between the hybrid and the factorized resiilts. 

We consider a process where an energetic projectile scatters off a static target. The wave function of the incoming 
projectile is an eigenstate of the QCD Hamiltonian. When calculated in the perturbation theory it can be represented 
as 



(2) 



where \v) is the zeroth order wave function (an eigenfunction of the free Hamiltonian), and is a unitary operator 
which diagonalizes the QCD Hamiltonian in perturbation theory 



The gluonic state immediately after scattering is 



(3) 



(4) 



where S is the eikonal scattering matrix for the projectile partons which propagate through the static target fields. 
The number of produced gluons is then given by 



dN 



1 



:{v\n^S''Q.a\k, k+)a{k, k+)n'' Sn\v) 



(Pkdk+ (27r)3 

Our first goal is to find the operator fl. We start with the light-cone Hamiltonian of QCD: 

H = l^^^^^d'z(^llL-{k+,z)IL-{-k+,z) + \G^^{k+,z)G:^i-k+,z)j 

where the electric and magnetic pieces have the form 

1 



(5) 



(6) 



G:^{x-,x) = d,A^{x-,x)-djAt{x-,x)-gr''^A^{x-,x)A%x-,x) 
Our convention for the covariant derivative is 

Jjab^b ^ f^g.^ab _ gf^cb^c^ ^b 



(7) 



(8) 



We arc; working in the light cone gauge hence A'^ = and as usual other light cone component of the vector potential 
A~ is expressed via the solution of Maxwell's equations as A^ = — The transverse components of the vector 
potential A^ which are the only dynamical degrees of freedom are expanded in the standard way in terms of the 
creation and annihilation operators 



Atix-,z) 



°° dk+ 1 



27r y/2k+ 

where the creation and annihilation operators satisiy the canonical commutation relations 

'atik+,x), y)l = 27r5«?'(fc+ - p+)S\x - y) 



(9) 



(10) 



We will calculate gluon production to the leading order in the coupling constant, and we therefore require to know 
the Hamiltonian only to first order in g. After some algebra we find 



H = Ho + Hi 



(11) 



Ho = [ ^af{k+,k)a'j{k+,k) 
Jk,k+>o 2k+ ^ 



k,p,k+,p+>o ^J2k+p+{k+ +P+) 



k+^' 



a!l{k+, k)a''j ip+ , p)af {k++p+,k+p) 



+ ■ 



-J^kja^{k+,k)a'i{p+,p)af{k++p+,k+p)'^+h.c. 
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where the integration measure is understood as 



dk+ 



and 



As f2 is a unitary operator, we define Hermitian operator G by 



-iG 



(12) 



The unitary operator $7 as discussed above is the operator that diagonalizes the Hamiltonian. To first order in the 
coupling constant the eigenvalues of the Hamiltonian arc those of the free Hamiltonian Hq. To this order we have 



Thus the operator G is determined from 



n'^HVl = H -i[H,G] =Ho 



i[Ho,G] = H, 



This immediately gives 



G = -gr'" [ . ^ = 

Jk,p,k+,p+>o \/2k+p+{k+ +P+ 



k+ 



) ^p+k ~ ~ l^k 



a\{k+,k)aUp+,p)af{k++p+,k+p) 



+ 



p+ + k+ 



kja1{k+,k)a^{p+,p)af{k++p+,k+p) \ + h.c. 



with 



(13) 
(14) 

(15) 



co{k) = 



2k+ 



(16) 



Gluon Production 



The number of produced gluons to leading order in the coupling is given by 



dN 



(Pkdk+ (27r)' 



S'fG- GS'' 



al\k+,k)alik+,k) 



GS-SG 



\v) 



(17) 



Here the factor j^jj is due to our normalization of the creation and annihilation operators eq.(lO). For simplicity 
of the calculation wc will assume that the longitudinal momentum of the observed gluon is (at least) slightly smaller 
that the momentum of gluons in the state \v), although in fact our formulae will be valid in a more general case. 

The calculation of the matrix element is straightforward. The S matrix operator acts as a color rotation on all 
gluon creation and annihilation operators in coordinate space 



(18) 



Since by assumption there are no gluons with longitudinal momentum fc+ in the state \v), one of the creation operators 



in the operator G in the amplitude 



GS-SG 



\v) must be at momentum fc"*" and is "contracted" with a{k~^) in the 

observable. This then leaves us with (apart from the various factors of S, and omitting for simplicity transverse 
dependences) expectation value of the type 



{v\a^p+ + k+)a{p+)a^q+)a{q+ + k+)\v) = 5{p+ - q+){v\a^p+ + k+)a{p+ + k+)\v) 

+ {v\a''{p+ + k+)a\q+)a{p+)a{q+ + k+)\v) 



(19) 



The second term involves a two particle density in the state \v). It is suppressed in the leading twist "'partonic"' 
approximation. Since we keep to this approximation in the present paper, we neglect this term. We note that in the 
soft approximation, where the gluon production is given by the kx factorized expression of [16] this term does indeed 
give a non vanishing contribution. We will make explicit connection with the soft approximation later. 
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Keeping only the first term in eq.(19) and reverting to coordinate space, where the S'-matrix is diagonal we obtain 
dN 



dPkdk+ 



(2^) 

^2 



1 



(^1 



1 



al\k+,z)al{k+,z) 



GS-SG 



^ik{z—z)-\-ipv-\-iqz—i{p-\-q)u—ipv—iqz-\-i{p-\-q)u 



(27r)3 7V2 _ 1 y fc+ 
2 



tr 



\v) 



(20) 



To arrive at this expression we assumed that the projectile state is color and rotationally invariant, so that 



Changing the integration variables 



1 



(21) 



p + q = K 
p + q = K 



(22) 



and integrating over w, cj, k, k, u and u one obtains 



dN 



1 



1 



d'^kdk+ 27r2 (27r)2 7V2 -ij^ ^ k+ 



^_}_^ik{z-z) 



(1-0 



(i-0' + e' + (i-0'^ 



{v - Z)i {v - Z)i 



X tr 
k+ 



{v — z)"^ (v — z)2 
^^((1 - Ov + iz)T-S{{l - + ^z) - SlT-S, 



,k+ 



,k+ 



(23) 



Throughout the rest of this paper we will continue using the notations u and u to make notations less cumbersome, 
however it should be understood that they are not independent variables, but rather as shorthand for 



u = {1 - + ^z; u={l-^)v + ^z 
To get some intuition about this expression we first consider the soft limit. 



(24) 



B. The soft limit 



The soft limit corresponds to the situation when the longitudinal momentum of the observed gluon is much smaller 
than the momentum of the gluons in the valence state \v). Taking ^ — )■ we obtain the this limit 



dN 



a., 1 



1 



d^kdk+ 7r2 (27r)2 N"^ ~ I J k+ {v - zY {v - zf 



(25) 



a. 1 



Nr. 



The kx factorized form reads ([16] and also [15], [17] corrected for typos) 



dN 



a.c 1 



d^kdv ~ 7r2 (27r)2 J ^ r^^' + ^''^^ ^'^^ ^"^M {z - v)^ {z - vf ^^"^^^ 



as 1 f 
(27r)2 j ' 



,ik(z — z) 



(26) 



where the last equality follows from color neutrality of the hadronic state. The color charge density operator here 

2-K 



Pv — I the leading twist approximation the correlator of the charge density operators 
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is local in the transverse space. Physically this is the case since in this (parton model) approximation there is only 
a small number of gluons in the hadron, and there are no correlations between different gluons. For a color singlet 
hadronic state we therefore have 

{p%v)p''{v)) = 6\v - v)N,{ J ^c4-{p+,v)anp+,v)) (27) 

Thus in the leading twist approximation eq.(25) is indeed equivalent to eq.(26). 

It is customary to define the transverse momentum dependent gluon distribution in terms of the gluon distribution 
function fg{x,Q) as 

(^'Q = ^73^) dH{at\p+,u)a1{p+,v)) = J dH J fle'''<^-^U{p,b;x) 

1 

« j d^b dp^(l){p,b;x) (28) 

where b = ^^^^ . In the soft limit the color charge correlation function and the scattering amplitude are then expressed 
in terms of the projectile and the target distributions as 

{p''{v)p'^{v)) = d^pe^^^-^^p'MPM (29) 

tr[l - S\v)S{v)] = 2TTasN^ I dPp^P-^-"^ \4,t{j),1)) (30) 

In terms of the transverse momentum distribution the single inclusive gluon spectrum in the soft limit has the familiar 
kt factorized form 



dN 
d'^kdrj 



S 



- 1 Ji 
a,N, 



I ■ k 



+ kf {i + kfk^ (i 



kf fc2 



,j>T{l + k)cj>p{l) 



(31) 



where we have assumed translational invariance in the transverse plane. Here S is the total transverse area of the 

collision and Y denotes the total rapidity difference between the projectile and the target in the process. 

In the limit of large momentum of the produced gluon fc 3> Qs, ^qcd the momentum integral in cq.(31) is dominated 
by two regions of momentum space. 

In the first region, I <^ k the dominant term is the first term in eq.(31) (which corresponds to the first term in 
eq.(25)). In this kinematics the incoming projectile gluon has a small transverse momentum (in accordance with the 
simple parton model picture) and it acquires a large transverse momentum due to elastic scattering from the target 
field. We will refer to this contribution as elastic: 



dN 
d^kdr] 



- elastic 



asNc 1 



SMI) 



(32) 



The final states that correspond to this contribution have a single high px gluon at forward rapidity. The balancing 
transverse momentum is carried by another gluon kicked out of the target by recoil, and it resides at a very different 
rapidity, close to the target. 

There is however another contribution which is equally important at the leading twist. This comes from the 
momentum range I = k + q with q <^ k. Changing variables from I to q the other contribution is clearly just the 
mirror image of eq.32). For reasons explained below we will refer to it as inelastic contribution 



dN 
d?kdr] 



inelastic 



asNg 1 

m-ik^' 



sMq) 



(33) 



q<Q~k 



In this kinematics all terms in eq.(31) are equally important. This contribution corresponds to a projectile gluon 
coming in with large transverse momentum in the wave function and subsequently scattering with small momentum 
transfer. The scattering practically does not add to gluons transverse momentum, but decoheres the gluon from the 
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incoming hadronic wave function. One can naturally ask, how do high transverse momentum gluons find themselves in 
the projectile wave function. The answer is, that they are always there as "'unresolved'" components of the '"parton 
model" ' gluons. A low pt gluon can split via a standard DGLAP evolution into a two gluon state with large relative 
transverse momentum. The " 'parton model" ' gluons are therefore not point like objects, but rather composites, which 
at first order in ag contain an admixture of a two gluon state. 

We stress that it is not the collinear part of the DGLAP kernel that is responsible for this structure. The coUinear 
emission contributes to multiplication of low momentum gluons in the wave function (" 'low" ' momentum here techni- 
cally means momentum lower than that imparted by the scattering). The splitting however contains also an ultraviolet 
contribution, which produces gluon pairs with large relative transverse momentum, which sit close to each other in 
the impact parameter plane. It is the presence of this compact two gluon configuration that makes a projectile gluon 
behave as a composite object. When such a composite parton scatters inelastically off a soft target field, its different 
components can be put on shell, emerging as high pt partons in the final state. Correspondingly, the structure of the 
final state incidentally is quite different from those that arise from the contribution eq.(32), as both high pT partons 
from the projectile wave function end up close to forward rapidity. The inelastic contribiition therefore takes into 
accoimt production of forward dijcts with large pt- 

To demonstrate that this is indeed the correct interpretation of eqs. (32,33) consider scattering of a "'composite"' 
projectile gluon state 



\in >= 



\v,a> +gT^J{v,z)\v,h; z,c> 



(34) 



The amplitude of finding a two gluon component in the soft limit should be identified with the Wezsacker- Williams 
field f{v,z) oc l^jEfyr- After propagation through the target 



lout >- 



S'^\v)\v,b> +gT,lfiv,z)S'"'{v)S-%z)\v,d; z,e> 



(35) 



Since the two gluon state admixture is perturbatively small, to leading order the elastically scattered state is simply 



\0Ut>elastic = S''\v) 



\v,b> +gT^J{v,z){z)\v,c; z,d> 



S''\v)\v,b> +gT,lS"'{v)S''%v)f{v,z)\v,c; z,e> 



(36) 



The amplitude of elastic scattering with momentum transfer k is then given by / (five'^''""[S'^''{v) — S""^]. We now 
square it, sum over the color indices of the final state, average over color in the initial state and multiply by the 
number of gluons in the incoming wave function with small transverse momentum. This gives the number of gluons in 
the final state produced via elastic scattering of gluons in the projectile, and is precisely the contribution of eq.(32). 
The inelastic component of the final state is 



\out >inelastic= \out > - \out >elastic= gThJiv, z)S {v) 



S'"'{z)\v,d; z,e > -S''''{v)\v,d; z,e> 



(37) 



We assume that the scattering matrix S{v) is a slowly varying function of the coordinate. Then the large transverse 

momentum of the gluon in the final state can only arise from a large relative momentum between the two gluons at 
V and z, that is from the high momentum component of the amplitude f{v, z). The amplitude of finding two gluons 
with large transverse momentum in the final state is given by 



9n 



d'^ze'''-^^ 



-v) 



f{v,z)S'>^iv) 



(38) 



Squaring this, summing (averaging) over the final (initial) color indices, multiplying by the density of gluons at v and 
integrating over the face of the projectile hadron we recover precisely the contribution of eq.(33) in the form of the 
second line of eq.(25). Thus we conclude indeed that the two leading twist contributions correspond to elastic and 
inelastic scattering of a projectile gluon on the target field. 

At high Pt, both contributions are of the same order of magnitude. The probability to find a low pr gluon in the 
projectile is of order unity, but then the probability of scattering with large momentum transfer is small, of order a^. 
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On the other hand, the probabihty to find a high pt parton in the incoming wave function is of order a^, however 
the probabihty of it scattering softly off the target is of order one. If one assumes that the projectile and target wave 

2 

functions have standard perturbative behavior, (j) = ii^, one finds 



dN 
cPkdrj 

dN 



aslJ.%HT In 



elastic 



d^kdrj 



QCD 

2 



2 2 1 

= asUpjjLrp In 



inelastic 



(39) 



where we have assumed that the perturbative behavior for the target kicks in at momenta above Qs . If the energy of 
the process is large enough so that the target distributions manifest extended geometric scaling[18], the Ql in the last 
equation should be substituted by a higher scale which marks the upper end of the geometric scaling window. At any 
rate, it is clear that at parametrically large transfer momentum the two contributions are comparable, and both must 
be kept. In the application we have in mind the transverse momentum is probably not much higher than Qs and so 
the inelastic contribiition does not have a logarithmic enhancement. However the logarithm in the elastic contribution 
is also not very large, perhaps a factor of 3 or 4. Thus whether the inelastic contribution can be neglected or not is a 
numerical question, and one would be well advised not to through it away prematurely. 



C. Back to the hybrid formaHsm. 



We now return to eq.(23). Our aim is to identify the two contributions described above in this more general formula, 
and to write the leading twist result in as simple form as possible. 

The elastic contribution as before corresponds to the region of the phase space where all the transverse momentum 
of the produced gluon originates from the momentum transfer from the target. This comes from the product of the 

two last terms in the brackets in eq.(23). Assuming that the momentum in the rest of the expression is small is 
equivalent to take z = z everywhere else apart from the scattering amplitude S'^{z)S{z) 



dN 



elastic 



asNc 1 
27r2 (27r)2 

-d. 



(i-0' + C' + (i-e)V 



eiK-s)tr[S\z)S{z 



1 



{v - z)2 27r^ 



(40) 



(27r)2 



where 



dp2 

p2<Q2 2p2 



g/g 



^Pg/giO^FfgC-^,p')NA{k) 



2N, 



(1 



0' + ^' + (i 



(41) 



is the standard gluon-gluon splitting function. As explained in [13] the first line is just the DGLAP contribution to the 
evolution of the gluon distribution. The " 'parton model" ' term is not present in our explicit formula because we have 
for simplicity assumed that the rapidity of the observed gluon is lower than that of any of the valence partons, and 
thus such gluons can be present in the wave function only as a result of evolution. This assumption is of course not 
necessary. Relaxing it restores the parton model contribution to the elastic gluon production. The elastic term then 
simply becomes the contribution discussed in [13] (we do not include yet the fragmentation function contribution) 



dN 
d'^kdrj 



elastic 



;XFfg{XF,Q^)NA{k) 



(42) 



with N^ik) = J d'^{z — z)e'^^'^^ ^Hr[S\z)S{z)]; and the factorization scale Q must be chosen so that it is of 

c 

order but smaller than the external momentum k. 

To extract the inelastic term in the leading twist approximation we note that it arises as the leading order expansion 
in powers of [a; — z| and \x — z\, since the separation between the two gluons in the wave function must be much 
smaller than the typical variation scale of the scattering amplitude S{x). Referring back to eq.(23) we write 



{V- 



- Sy 
Sy - 



- ^{v - z)idiSy 
{v - z)jdjSy 



(43) 
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Then for the ampUtude in eq.(23) we find 



tr 



»3i7.-i i~>U ~ '-',1-1 i~>Z 



'^z '^'^ 



= N,iv-z),iv-z),{il-0^+e}tr 
= ^{v-z)-{v-z){il-0'+e}tr 



di dj Sjj 



-2(,{l-i){v--zUv-z),tr 
- C(l - -z)-{v- z)tr 



T'^S.diSlT-S.djSl 



T-S^diSlT-SydiSl 



(44) 



where the last line strictly speaking holds only after averaging over rotationally invariant target for which {diS^djS) = 
^Sij{dkS^dkS). One can further simplify this expression noting that SydiSl = ■^T'^tr[SydiSlT°'] i.e. it is a '"pure 
gauge" ' vector potential. Then using simple color algebra we find 



tr[SdiS''T"SdiS''T"] = -^trldiS'fdiS] 



(45) 



So that we can write 



tr- 



N. 



SiT-Su - SlT-S, 



N^{l-^ + e}iv-z)-iv-z)tr 



(46) 



where we have given the final answer in terms of the fundamental representation matrices Sp- It is clear from our 
derivation that in the above expressions the target field average should be understood as calculated with resolution Q, 
where just as in the clastic piece, the factorization scale Q is of order of, but smaller than the transverse momentum 
of the observed gluon Q < k. 



tr 



diSF{v)diSUv) 



tr 



(47) 



The leading twist part of the inelastic contribution can therefore be written as 
dN 



dPkdrj 



a. 1 



1 



inelastic 



27r2 (27r)2 iV2 _ i 



'^-'^ ^-%^-'^ (v-r)-iv-z)tr 

{v — zy [v — zy 



diSF{v)d,sy{v) 



X Ni 



(1-0 



(48) 



This expression can be rewritten in terms of the gluon distribution. To leading twist the dependence of the vacuum 
average (a^^(^, M)a^(^, u)) can be substituted by {^j v)C'']{^tv))q. Performing the Fourier transform we then 
obtain 



dN 
d?kdr} 



Ois 1 d^ 



inelastic 



d^p 
(27r)2 



p'NFip) (49) 



Of course, in order to calculate the spectrum of produced hadrons we have to include the gluon fragmentation 
functions. We assume as always, that produced gluons fragment into hadrons independently Taking this into account 
our result for production is 



dN 
d'^kdrj 



^ dz 
—Dh/g{z,Q) 



xifgixi,Q^)NAix2,-,b = 0) 



(50) 



a, N;^ z 



_fp_ 

{2nf 



p^NF{x2,p,b = 0) 
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where 



NF\k,b 



2 

z + z 



(51) 



^ J (f{z-z)e"'<^-'hr[SU^)SF{z)] 



and the longitudinal momentum fractions (neglecting the hadron mass) are 



Xf 



k 



Xf 



, e'; x-L 

\/snn z 



X2 = x\e 



-2r) 



(52) 



Eq.(50) is our result for hadron production in a toy theory that does not contain quarks. This is obviously not a 
good approximation to reality especially at forward rapidities, where the quark contribution must be the leading one. 
We now turn to generalizing the previous discussion by including the quark contribution. 



III. TAKING CARE OF QUARKS 

We now include the quark piece into the light cone QCD Hamiltonian 



H„ = J+- 



where 



and as usual 



P'' = -iDi"; Di=di+ igr^AI; Ifr = ^'Pi 



Here t" are the generators of SU{N) group in the fundamental representation. 

The dynamical spinors ^1/+ are defined in terms of creation and annihilation operators as[19] 

^l{x-,x) = [ (fkdk+ «;,6^(A;+,fc)e-*'=''^"^«^+«;_«4«(fc+,fc)e**^''^"+''=^ 

(27r)'^ Jk+>o L 

Here bg {dg) are the quark (antiquark) annihilation operators, and the basis spinors Wg are 



Wi 



, W_i 






\1/ 



(53) 

(54) 
(55) 

(56) 



(57) 



The quark algebra is 

{b,ik+, k),bl,{k'+, k')} = {4(fc+, fc),4(fc'+, fc')} = (2^)3<5(fc+ - fc'+)(52(fc - k')S,s' (58) 
When expanded in terms of quark, anti-quark and gluon creation/annihilation operators Hq can be written as 

Hg = Ho+ Hi (59) 

where the free Hamiltonian Hq is 



Hn 



(27r)3 



d^kdk^ 



k+>0 



2k+ 



b'^\k+,k)b':{k+,k) + df{k+,k)d'^{k+, k) 



(60) 
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and 0{g) Hamiltonian Hi is 



f 2p+ + 



Pt Qt 
p+ q+ 



p+ q+ 



k+>o \/2q+ { P+q+ 

bl?{p++q+,P + q)b^s{P^,PXjiQ^,q) + d'Xp+,p)dl^ip++q+,P + q)a'^,iq+,q) 
1 



+ h.c. 



q+ -p+ 



_p+ +q+ 

X g)6f (p+,p)at'^(p+ + + + h.c. 



Pi_ _ (h_ 

p+ q+ 



Pt_ _ <h_ 

P+ q+ 



Following the same steps as before, we find the operator that diagonahzes the Hamiltonian to first order 

O = e*^ = 1 + iG - + ... 



with 



G 



'p,q 

+ 



1 



2W2q+C0p+q 



_J r _ . 2p+ + q+ 

- Wp - \ p+ + q+ 



Pi Qi 
p+ q+ 



diid..'.x - e*Vf- — — 



p+ q+ 



bl?{p++q^,P+q)bUp^,p)<{q+,q) + d':,ip+,p)dl^{p++q+,p + q)a'^Jq+,q) 



1 



1 



■ q -p 

1 

2y/2{q+ +P+) i^p+q - UJp - OJq [ p+ + q+ 

X 9T^pd'^,iq+,q)bl{p+,p)a]'{p++q+,p + q) 
+ h.c. 



Pi qi 

p+ q+ 



Pi_ _ qi_ 
p+ q+ 



Again, following the same steps as before we calculate the number of produced gluons 

= e'''^^-'\v\[S^G-GS^]4''{k+,z)at{k+,z)[GS-SG]\v) 

{v - z)i [v - z). 



d'^kdk+ 



as 1 
27r2 (27r)2 



k+ [v — zy [V — zY 



- S^-^z) 



- S^^^iz) 



T{^,u)K{'^,u)) + {dr{^,u)dU'^,u)) 



(61) 



(62) 



(63) 



(64) 



In the leading twist approximation there are, as before, contributions corresponding to elastic and inelastic partonic 
scattering. The elastic contribution can be written as 



dNg 
d^kdri 
k+_ 

277^ 

a 



'^)^ Jxf C 



a 

elastic ^TT^ (2„; j^^ 

k+ ..^.k+ ,....A;+ 



{bl^{^,u)K{^,u)) + {dl^{'^,u)d^,{^,u)) 



o^s I [ dp f' 



{V - Zf {V - 2)2 Nl - 1 



NA{k) 



S\{z)Sa{z) 



(65) 



with 



P9I9- 27V, 



iV2 - 1 1 + (1 - if 



(66) 
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This takes into account the DGLAP correction to gluon distribution due to quark(and antiquark)-gluon spUtting. 

When written in terms of gluon distribution, this term is aheady contained in eq.(42). The inelastic contribution is 
given by the region of phase space where the momentum transfer from the target is small. Hence one can expand the 
eikonal factors S{u) around S{v). The result is 



dN„ 



cPkdr] 

k+ 



27r^ 



inelasUc ^^'^ {^nf ^ ^ ^ ^ _ ,-)2(^ _ ,)2 

-tr djSl{v)djSF{v) 



2 1 



(67) 



iV2 1 



27r2 m - 1 fc4 



1 + 



/,(y,Q) + /,-(y,Q) 



{2nf 



Similarly one can also calculate the number of produced quarks 



dNn 



1 



d'^kdk+ 



(27r)3 
a,, 1 



Jk{z-z), 



;| [SG - GS''] h\''{k+,z)h'l{k+, z) [GS'^ - SG] \v) 



+ 



27r2 (2 

as 1 
27r2 (27r)2 



1 / e*M-?) 1 (1 + ^2 



{v - z)-{v- z) ,^^gjk+ -^,„(k+ . 
(7; - z)^(i) — z)^ ^ ^ ^ ^ 



Sf{z)t^ 



ik{z — z 



^ (f^ + (1 - ;::|;::$ (»r(f.^)»;(f.»)> 



X tr 



Spiuys'j.iu) - spizysUv) 



(68) 



The elastic contribution is 



dNg 
d?kdr} 



as 1 



elastic 



47r2 Nc (2 



k+ 



k+ 



,k+ 



1-^ (w-^)2(t;-z)2 7V^ 



tr{sy{z)SF{z))—M''{^,u)b'l{—,u)) 



as 1 /■ ^ 



(69) 



with 



^l±^- P, =i(.2 + (l_.)2) 



(70) 
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and the inelastic term is 



(Pkdrj 



asN^ 1 



inelastic 



87r2 (2 



(i + a-(i-OV 



[{v-z).{v-z)r 1 



1,+ 1.+ 1.+ 



167r2 (27r 



1 / |e^M-.)(^2^(l_^)2) 

7r)2 7^;^ C 



[(l_^)2+^2]_2C(l-€) 



(I) — zy{v — 2;)2 



as I f d^P 2at r \ '^^ 



(71) 



iV2 



7V2-1 



p2<Q 

1+e- 



2 



iV2 



^ iV2-l 



For antiquarks the calculation is identical. The elastic contribution is 



dNg 
dPkdr] 

J elastic 

and the inelastic term is 



Ola 1 





di 









^f(fc) (72) 



dPkdr] 



inelastic 



2 (1-C)^ 



(73) 



7V2 
7V2- 1 



1 + ^^ 



P,/,(C)/«-(y,Q) + ^ 



It is convenient to introduce the following functions, which we will dub "inelastic weights" 



m-1 



7V2 



iV2 

1 



2 



Wq/giO = Wq/giO 



iV2- 1 



Pq/giOfgiY'Q^ 



(74) 
(75) 
(76) 
(77) 



In terms of these functions one can write the inelastic contribution to production of i'th partonic flavor as 



dNj 
d'^kdr] 



- inelastic 



Finally, including the effects of parton fragmentation, the result for hadron production is 

dNu ^ dz k k 

' x,f,ixi,Q^)NA{x2,-)Dh/g{z,Q) + ^gX,fgixuQ^)Np{x2,-)D^/g{z,Q) 



d^kdv (27r)2 J^^ z2 
dz as -2* 



z 

Xi 



iS^hl^^Q. (If^^'^^^^'^^)^^/ j^j=,,^,gWi/AOPi/AOfj{j,Q)DH/q{z,Q) (79) 

where the momentum fractions Xi and X2 are defined in eq.(52) and the inelastic weights Wi in eq.(74). This is our 
final result. 
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IV. DISCUSSION 



In this paper we have derived the complete leading twist expression for inclusive hadron production in the hybrid 
formalism. We have shown that in addition to elastic scattering terms first derived in[13], there are also terms that 
correspond to inelastic scattering of the projectile partons on low momentum components of the target field. These 
terms are given by the second line in eq.(79). We note that although the inelastic piece has an explicit factor of 
while the elastic contribution does not, the two terms at high kr are in fact of the same order in as- The reason is 
that at momenta k ^ Qs the dipole scattering amplitude Nj^i^p-f{k), which enters the elastic scattering term is itself 
of order as , while the integral of the amplitude appearing in the inelastic term is of order unity. 

The final states that correspond to the inelastic process are dihadron pairs where both hadrons are emitted at 
forward rapidity and have strong back to back correlation. Since both produced hadrons have large rapidity, such 
pairs with large transverse momentum are kinematically allowed only at large collision energy. Thus one might expect 
this contribution not to be of great importance in RHIC kinematics, however it may be sizable at LHC. 

In this context we believe that including this contribution in calculation a la [7] should produce faster approach 
of nuclear modification factor RpA to unity at large transverse momenta. Here we wish to elaborate on possible 
role of saturation in the resiilts of [7]. As we have noted above, as long as the transverse momentiim is above Qs, 
saturation should mainly affect the inelastic production piece. This contribution involves the target distribution 
/p2<-Q2 i^J')-i p'^Nf{P: X2) oc ftargetiQ,X2) and is thus directly sensitive to saturation effects which suppress the con- 
tribution of small momenta p < Qs to the integral. The elastic production probability (first line of eq.(79)) depends 
only on N{k) at large momentum. Naively one expects that this part is unaffected by saturation in the evolution. 
This does not necessarily mean that the RpA calculated using only this contribution (as done in [7]) should be equal 
to unity, but rather that any departure from unity is the effect of a non scaling initial condition. This should be 
true if the transverse momentum of the measured particle is in the so called '"geometric scaling'" window, where the 
anomalous dimension is finite, since geometric scaling is not a result of saturation physics but rather of the linear 
BFKL evolution of the gluon density [18]. To be a little more precise, recall that solution of the BFKL equation above 
the saturation scale has the form 4>BFKL{k, Y) oc [Q s{Y) / k]'^~^^ where 7 is the anomalous dimension. The anomalous 
dimension is a slowly varying function of transverse momentum. It is almost constant in a wide window of momenta 
above Qs, but nevertheless vanishes asymptotically as fc — > 00. It also weakly depends on rapidity Y . The saturation 
momentum Qs is defined withing the BFKL solution per se as the momentum at which the scattering amplitude is 
of order one. Within leading order BFKL solution Qs(Y) = Qoexp{XY} where Qq is the soft nonperturbative scale 
which characterizes the initial condition (j)o{k,Y = 0). In the case of calculations of [7] this would be the initial 
saturation scales, Qop for the proton target and Qqa for the nuclear target. The nuclear modification factor RpA 



within a BFKL calculation would then be 



QsAjY) 

QspiY) 



2-27(F) 



QoA 



2j{Y) 



(80) 



with the identification Ncoii — Qqa/Qop- Within the running coupling calculation the saturation scale is not a simple 
exponential of rapidity and thus the explicit expression for the nuclear modification factor and the rapidity dependence 
is somewhat different. It nevertheless remains the case that as long as the initial conditions for proton and nucleus do 
not simply scale with A^/^ at all momenta, 0p(fc, Y = 0) ^ A^^^(j)A{k, Y = 0), the linear BFKL evolution produces a 
nuclear modification factor different from unity and slowly varying with rapidity. It is an interesting question whether 
the numerical results of [7] are consistent with BFKL, or whether saturation effects in the evolution nevertheless give 
a significant contribution to RpA- 

Finally we note that the final states that contribute to the inelastic production arc precisely the states which give 
the bulk of the contribution to the dihadron correlation function considered in [10]. The calculation of [10] does not 
address the estimate of large uncorrelated background of produced particles. The small " 'signal to background" ' ratio 
is indeed a very pronounced feature of the data [9], [20]. In this respect it would be interesting to calculate both in a 
unified framework discussed here. We note that our earlier discussion suggests that the saturation has two distinct 
effects on the correlated dihadron production. First, as discussed in [17], [21] and [7], the back-to-back correlation is 
weakened due to independent momentum transfer from the target to each one of the produced hadrons. This does 
not reduce the number of hadrons produced at forward rapidity, but reduces the correlation between the direction of 
their transverse momenta. Another distinct effect is that the dihadron production probability is suppressed by the 
effect of saturation on /p2<Q2 ■^:j^P^Nf{p,X2), thus reducing the ratio of the correlated signal to the total number 
of produced particles. 
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